Multitype branching processes with immigration 
in random environment and polling systems 



Vatutin Vladimir* 
Steklov Mathematical Institute, 
Gubkin street 8, 119991, 
Moscow, Russia 
E-mail: vatutin@mi . ras . ru 



Abstract 

For multitype branching processes with immigration evolving in a ran- 
dom environment and producing a final product we find the tail distri- 
bution of the size of the final product accumulated in the system for a 
life period. Using this result we investigate the tail distribution of the 
busy periods of the branching type polling systems with random service 
disciplines and random positive switch-over times. 



1 Introduction 

Branching processes with and without immigration arc powerful tools in study- 
ing various models of queueing systems (see, for instance, [25], [26], [27], [28], [35], 
[36], [41] and [43]). In this paper we use multitype branching processes with 
accumulation of a final product and immigration which evolve in random envi- 
ronment (MBPFPIRE) to study the tail distribution of busy periods of a class of 
polling systems in which input parameters, service disciplines and the distribu- 
tions of switch-over times vary in a random manner. This article complements 
the results of paper [16] established for polling systems with zero switch-over 
times. 

The paper is organized as follows. In Sections 2-6 we established various 
results related to MBPFPIRE. These results are applied in Section 7 to inves- 
tigate the tail distribution and moments of the branching type polling systems 
with positive switch-over times. 
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2 Branching processes in random environment 
with final product 

Let s : = (si, . . . , s m ) € [0, l] m be a m-dimensional variable, 



s k := s* 1 



fci€ N = {0,1,2,...}, 



and (£; <p) be a (m + l)-dimensional vector, where the vector £ := . . . , £ m ) 
has integer- valued nonnegative random variables as components and tp is a non- 
negative random variable, and let 



F(s;A) := E [s«e~ Av ] - 



E 



Jl J2 

6 1 *2 



s g [0, l] Tn , A > 0, 



be the respective mixed probability generating function (m.p.g.f.). Denote F : 
{F(s;A)} the set of all such m.p.g.f.'s and let 

F m := {F(s;A) = (FW(s;A), . . . ,F( m >(s;A)) : F«(s;A) 6 J, t = 1, ...,m} 

be the m-times direct product of F. Let, further, 



s^ 1 <j£ 2 . . . <j£« 



F(s;A) e F} 



F := {/(s) = F(s;0) = E 

be the set of all ordinary probability generating functions (p.g.f.'s) and 

FJ l := {f(s) = (/W( S ), . . . , /W( S )) : /W(s) 6 Fo, i = 1, ...,m} 

be the set of all m-dimensional (vector- valued) p.g.f.'s. 

Let, further, (to ip) be a (to + l)-dimcnsional tuple where the vector ?/ := 
{rji , . . . , ?7 rn ) has integer- valued nonnegative random components and ip is a non- 
negative random variable and let 

G(s; A) = E [ s"e- A ^] = E [ sf sf ■ ■ ■ e"^] , s e [0, 1]™ , A > 0, 

be the respective m.p.g.f. Denote Q := {G(s;A)} the set of all such m.p.g.f.'s 
(which, of course, is equivalent to F). Let, further, 



g o :={g(fi) = G{s;0) = E[sf 



3 m J ' 



G(s;A) e g} 



be the set of all ordinary probability generating functions (p.g.f.'s) . 

Assume that a probability measure Q is specified on the natural cr-algcbra A 
generated by the subsets of J rm xg. Let H(s;A), H (s;A), Hi (s; A), ... , Hfc(s;A), . . . 
where 

H(s;A) := (F(s;A); G(s;A)) = (f«(s;A), . . . , F<"0(s;A); G(s;A)) , 

and 



H„(s;A) := (F„(s;A); G„ +1 (s;A)) = ( F« (s;A), . . . , F^(s;A); G„ +1 (s;A) ) , 
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n = 0,1,2,... be a sequence of vector-valued m.p.g.f.'s selected from T m x Q 
in an iid manner in accordance with measure Q. The sequence {H ra ,n G No} 
is called a random environment. With each m.p.g.f. Fn^ (s;A) we associate a 

(&i( n ),&2(n), . . . ,iim(n)) and a random 



random vector of offsprings £j(n) := 
variable ifi(n) such that 



^ 4) (s;A) 



E 



8 £<(n) e - 



Ac^i (n) 



E 



and with each m.p.g.f. G„(s;A) we associate a random vector of immigrants 
77(71) := (771(72), rfein), . . . , r/ m (n)) and a random variable V>(n) such that 



G„(s; A) := E [s^V^M = E 



7)1 (n) 7)2 (n) 



sl 2(n) e" 



Ai/>(r 



= G(s;A). 



Now we may give an informal description of a multitype branching processes 
with accumulation of a final product and immigration which evolve in random 
environment (MBPFPIRE) 

T(n) = (V(n);6(n)), neN , 

which may be treated as a process describing the evolution of a m— type popu- 
lation of particles with immigration and accumulation of a final product. 

Given an environment {H„,ti e No} the starting conditions for the process 
are: a vector of particles V(0) = (Vi(0), . . . , V^,(0)) (may be random or equal 
zero), where Vi(0) denotes the number of particles of type i € {1, . . . , m} in the 
process at moment 0, and an amount 0(0) (may be random or equal zero) of a 
final product. All the initial particles have the unit life length and just before the 
death produce children and some amount of the final product independently of 
each other. A particle, say, of type i, produces at the end of her life particles of 
different types and adds some amount of the final product to the existing amount 
of the final product in accordance with m.p.g.f. (s;A). In addition, a ran- 
dom tuple of immigrants specified by the vector 77(1) =: (771(1), 772(1), 77 m (l)) 
arrives to the system at moment 1, where 77^ ( 1 ) is the number of type j particles 
immigrating in the first generation of the population, and the final product of 
size t/>(1) is added to the system with the joint distribution specified by the 
m.p.g.f. Gi(s;A). Thus, the amount of the final product 0(1) accumulated in 
the system to this moment is 

rn Vi(0) 

6(1) = 6(0) + V(1)+^ <Pi(0,k), 

7=1 k=l 

where ipi(0, k) - is the amount of the final product added to the system at the 
death moment of the fc-particle of type i among those which existed at moment 
0. 

The newborn particles and immigrants entering the system at moment n = 1 
constitute the first generation of the MBPIFPRE, have the unit life-length and 
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dying produce, independently of each other and of the prehistory of the pro- 
cess, offsprings and final product in accordance with their types and subject 
to the m.p.g.f. F^ l \s;X),i = 1,2, ... ,m. In addition, at moment n = 2 immi- 
grants and some amount of final product described by a tuple (77(2); tp(2)) ~: 
(771(2), 772(2), n m (2); ifi(2)) is contribute to the process, which is specified by 
the m.p.g.f. G2(s;A). And so on. 

Note that in our model "0( n ) m ay be positive even if rj(n) = 0. 

Definition 1 A m-type Galton-Watson branching process 

T^(n) =T(n) := (V(n);0(n)) = (V^n), V m {n); 6(n)) , n 6 N 

with immigration and final product (92, VO in a fixed (but picked at random) 
environment {H„, n € No} is a time-inhomogeneous Markov process with the 
state space 

N " x R + := {z = (zi,...,z m ;w), Zi € No; w £ [0,oo)} 

defined as 

T(O) = (V(O);e(O)) = (v;0), 
E [ s v(n+i) e -Ae(„+i) i Ho; Hn . T(Q)j ; T(n) j = e -X0(n) (F „( s; A)) V(n) G„ +1 (s;A). 

Thus, 

T(n + 1) = (V(n + l);6(n + l)) 

m Vi(n) 

= (0; 0(n)) + (r?(n + 1); i>(n + 1)) + J2 J2 *)i ^ ' 

where (£i(n; fc); <^i(n; fc)) is a random vector representing the offspring vector and 
the size of the final product contributed to the process at the death moment of 
the k— th particle of type i of the n— th generation. Given the environment and 
V(n) the tuple 

(&(n;fc);<^(n;fc)), k = 1, 2, . . . , V l {n), i £ {l,...,m}, (r?(n+l), »» £ N ; 

consists of independent vectors and, moreover, for each n £ No and i £ {1, . . . , m} 
the vectors 

(&(n; fc); <^(n; fc)) , fc = 1, 2, . . . , Vj(ra) 

are identically distributed: (£i(n; fc); <^i(n; fc)) = (&(t5.); ¥'i( n ))- 
In the sequel we write for brevity for fixed n £ No 

Ph„ (•) := Ph„ (• I H„) , E H „ [•] := E [• | H„] , 

Ph (•) := Ph (• I H , H„) , E H [•] := E [• | H , H„] 
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or (where it will cause no confusion) 



Ph (■) := Ph (• I H , H n , ...) , E H [•] := E [■ | H , H„, ...] 

Similar meaning will have the notation Pf„ , Ep„ , Pp, Ef and so on. 

Observe that if 9(0) = and (p z (n; k) = 1, ip(n) = 0, n <E N , then Q(N) - 1 
is the total number of particles born and immigrating in the process within 
generations 0, 1, N — 1; if 

(Pi(n; k) = I | ^ ey (n; k) > t J , j>(n + 1) = 0, n € N 

for some positive integer t (here and in what follows I {A} stands for the indi- 
cator of the event A) and 9(0) = 0, then Q(N) is the total number of particles 
of all types in generations 0, 1, ...,N — 1 each of which has at least t children, 
and so on. 

Letting A = in Definition 1 we arrive to the definition of the multitype 
branching process with immigration evolving in random environment (MBPIRE) 
which we call the underlying MBPIRE for the initial MBPIFPRE. 

Definition 2 A m-type Galton-Watson branching process with immigration V(n) 
in a fixed (but picked at random) environment {(f n , flVi+i)) n G No} is a time- 
inhomogeneous Markov chain with the state space 

^ :={z = (z 1 ,...,z m ),z i eN } 

defined as V(0) = z and for n £ No 

E [s v (' l+1 )| f , 5l ,...,f„, ffn+1 ;V(0),...,V(n)] = (f„ (s)) V(n) g n+1 (s). (1) 

Thus, 

m Vi(n) 

V(n + l) = J2J2 £<(»;*) +*7(n+l). 

i=l fc=l 

where, given the environment and fixed V(n), the tuple 

£i(n;k),k = 1, 2, ...Vi(n); i = 1,2,..., to; rj(n + 1) 

consists of independent random variables. 

In the sequel we need the definition of the ordinary multitype branching 
process with final product which evolves in random environment (MBPFPRE) . 

Definition 3 A m-type Galton-Watson branching process 

R v (n) = R(n) := (Z(n); $(n)) = (Z^n), Z m (n); <P(n)) ,neN , 



5 



with final product tp in a fixed (but picked at random) environment {F„, n £ No} 
is a time-inhomogeneous Markov process with the state space 

N™ x R := {z = (z 1 ,...,z m ;w), z t e N ;w £ [O,oo)} 

defined as 



E 



R(0) = (Z(0);*(0)) = (z;vo), 

8 Z(n +l ) e -A*(»+l) | Poj . . . } Pn ; R(0 ), . . . , R(„)l = e -A*(n) (F^A)) 2 ^ 



Note that the initial value R(0) may be random. Thus, 

<£(n+i) = m+J2Y,I2^ l -> k ) ( 2 ) 

1=0 i=l fe=l 

where (£i(n; fc); tpi(n; k)) is a random vector representing the offspring vector 
and the size of the final product of the k— th particle of type i of the n— th 
generation of the process. 

Finally, excluding immigration and final product we obtain the definition of 
the ordinary underlying multitype branching process in random environment. 

Definition 4 A m-type Galton-Watson process 

Z(n) = {Z x {n),...,Z m {n)),n £ N 

in a fixed (but selected at random) environment {f„, n £ No} is a time-inhomogeneous 
Markov chain with the state space 

Nq™ :={z= ( Zl ,...,z m ), Zi eN } 

defined as 

Z(0)=z, E[s z (" +1 >|fo,...,f„;Z(0),...,ZH] = (f„ (s)) z(n) . (3) 

It is necessary to mention that the one-dimensional BPRE without immi- 
gration was introduced by Smith and Wilkinson in [39] and, in a more general 
setting in [G] and [7] and has been investigated by many authors (see survey [17] 
for a list of references up to 1985 and [1],[2],[3], [18]-[21],[23],[24], and [33] for 
some more recent results). 

The ordinary single- type Galton-Watson branching processes with final prod- 
uct where investigated by Sevastyanov [38] for the particular case tp(n; k) = 1, 
and by Grishechkin [27] (for the general ip(n;k)). MBPRE were analyzed in 
[ ] and [40] and MBPIRE were considered by [32] and [37]. MBPFPRE were 
introduced in [16] to study the tail distribution of busy periods of the branch- 
ing type polling systems with zero switch-over times. In the present paper we 
analyze MBPIFPRE and apply the obtained results to investigate properties of 
busy periods of the branching type polling systems with random and, generally 
speaking, positive switch-over times. 



G 



Let ej := (0, . . . , 1, . . . , 0) be a to— dimensional vector with zero components 
except the i-th equal to 1, and 1 are to— dimensional vectors all whose com- 
ponents are equal to and 1, respectively. 

Denote 

A I I \\ m I ®fn ( S ) 

A n = (a 4i (n)) ij=1 := 



ds o 



s=l 



(4) 



the mean matrix of the vector-valued p.g.f. f n , 

C„ := (E Fn ^ 1 (n),...,E F ^ 1 (n))' = ( ^}*' X) 

'(5) 

the mean vector of the amount of the final product of particles of the n-th 
generation and 



A=0 



dFt\l,X) 



A=0 



B„ := (Bi(n),...,B m (n))' 



dg n (s) 



dg n (s) 



D n := E G „V(n) = 



dsi 
dG n {l,X) 



ds ri 



dX 



A=0 



the respective characteristics for immigrants. 

By our assumptions the tuples (A n , B„ + i, C n , D n+ \) , neffo are iid: 

(A,, B n+ i, C n ,D n+1 ) = (A, B, C, D). 
For vectors u = {u\, . . . , u m ) and v ={v\, . . . , w m )' £ K m denote 



(6) 
(7) 



(u, v) := ^2 u i v -i 



fc=l 



their inner product. 

Foramxm matrix ^4 = (tiy) i J=1 and a to— dimensional vector u =(iti, 
introduce the norms 



and 



I4| 2 := 



\ 



Let, further, 



n, 



|u|| :=X>I 
i=l 



\ *=1 



1< Z < n, 



with the agreement that II n ,„ := E is the unit to x to matrix. 
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3 Statement of main results 



Denote T> := {x > : E \\A\\ X < 00} and, for a given x > 0, set 

s(x):= lim (E|K_ 1 ..-A |r) 1/ "= lim (E ||n ,„ \\ x ) 1/n (8) 

and let 

a:=s'(Q)= lim ^Elog \\A n ^ ■ ■ ■ A \\ = lim -Elog ||n , n || (9) 

n— too n n— >oo n 

be the top Lyapunov exponent for this sequence of matrices. 

It is known that the limits in (8) and (9) exist and, moreover, s(x) is a 
log-convex continuous function in T> (see, for instance, [34]). Put 

K := inf {x > : s(x) > 1} (10) 

and k = 00 if s(x) < 1 for all x > 0. Observe that s(0) = 1 and, therefore, 
k = if a > and « e (0, 00] if a < 0. 

In the last case (which will be our main concern) the series X)^Lo ^ ||Ho,r» ll 35 
converges if < a; < k and diverges if x > k. 

In what follows we call the underlying MBPRE subcritical if a < and su- 
percritical if a > 0. The same terminology wc keep for the respective MBPIRE. 

Now we formulate important statements concerning properties of MBPRE 
and MBPIRE. 

Let 

ft(f) := lim P f (||Z(»)|| =0|Z(0)=ei) 

n— > 00 

be the extinction probability of a MBPRE initiated at time by a single indi- 
vidual of type i and 

q(f) := ( qi (f),...,q m (f)) ■ 
Theorem 5 ([40]) If the mean matrices of a MBPRE meets the condition 

Elog+ \\A\\ < 00 (11) 
anc? there exists a positive integer L such that 



min (Al^Al-2 ■ ■ ■ A ) H > = 1 

l<2j<m J J 

and 1 < I < m sitcft 

E|log(l-P f (Z I (L) = 0|Z(0) = e i ))| < 00, 

tten, /or a specified by (9) 

1) a < impZies P (q(f) = 1) = 1; 

2) a > implies P (q(f) < 1) = 1 and 

P f ( lim n- 1 log ||Z(n)|| = a \ Z(0) = e. t ) = 1 - $ (f) (12) 
with probability 1 for 1 < i < m. 



The next statement deals with MBPIRE. 
Theorem 6 ([32] and [37]) Let a MBPIRE satisfy condition (11), a < and 

Elog+ ||B|| < oo. 
Then, for any v £ N™ the limit 

lim P (V(n) = v) =: £>(v) 

n— >oo 

exists and defines a probability distribution on W n . //, in addition, there exists 
L > 1 sucfi. that 

P(P g (V(L) = 0|V(0) = 1) > 0) > 
i/ien -D(O) > and, therefore, 

P(V(n) = i.o. ) = 1. 

Introduce the set 

U+ = {u = ( Ul , . . . ,u m ) £ R m : m > 0, 1 < i < to, ||u|| a = 1} 
and associate with the tuple (A n , C„ ), n £ No of iid pairs the series of vectors 

oo oo 

Ei := AiAi+x ■ ■ ■ A k ^C k = ^ U^ k C k , I £ N ; S := S . (13) 

fc=i fc=; 

Our main results are established under the following hypothesis. 
Condition T. There exist positive constants k and Kq and a continuous 
strictly positive function l(u) on U+ such that for all u £?7+ 

lim y»P«u,S) >y)=K l(u). (14) 

In Section 4 we list sufficient conditions imposed on the distributions of the 
pairs (A n , C„) which provide the validity of Condition T. These conditions are 
extracted from paper [30] where the behavior of the tail distribution of sums 
and products of random matrices were investigated. 

Let 

oo m Zi(n) 

$ := lim $(n) = $(0) + VVV ipAn\ k) (15) 

n— foo ^ — ^ — J — ^ 

n=0 i=l fc=l 

be the total size of the final product produced by the particles of the MBPFPRE 
up to the extinction moment (if any). 

The following theorem was proved in [16]. 

Theorem 7 Let a MBPFPRE satisfy the following hypotheses: 

1) the underlying MBPRE is subcritical and meets the conditions of Theo- 
rem 5; 



2) for n specified by (10) the following assumptions fulfill: 
if K > 1 then 




if k < 1 i/ien 



max E 

l<i<m 




(17) 



^ there exists S > suc/i i/iai < E<^ + < oo, i = 1, . . . , m; 

5) i/ie mean matrix (4) and the vector (5) are such that Condition T is valid. 

Then, as y — > oo 



The theorem has been used in [16] to investigate the tail distributions of the 
busy periods and some other characteristics for a wide class of polling systems 
with zero switch-over time. In the present paper we consider MBPIFPRE and 
apply in Section 7 the respective results to deduce similar statements for busy 
periods and certain other characteristics of polling systems with positive (pos- 
sibly random) switch-over times. To this aim we introduce the notion of life 
periods of a MBPIFPRE. 

Definition 8 We say that a branching process T(n),n £ No with m types of 
particles, immigration and final product has a life period of length T n initiated 
at moment N > 1 if (V(N - 1); Q(N - 1)) = (0; 6), \\rj(N)\\ > 0, and 



In what follows we consider without loss of generality the distributions of 
various characteristics related with the life period of length T := Ti, and study 
the tail distribution of the quantity 



- the total amount of the final product accumulated in the MBPIFPRE during 
the life-period which starts by 77(1 ) = z / particles at moment N = 1. Since 
our proofs (but not constants!) do not depend on the particular value of 77(1) 
we do not specify it explicitly in the subsequent arguments. 

Note that if 9 = and (fii(n; k) = 1, ip(n) = for all n G No then is the 
total number of individuals existing in the MBPIFPRE during the respective 
life-period. 

Now we are ready to formulate the main result of the present paper. 



P($>2/|Z(0) = z)~C(z)y 



c(z) e (0,00). 



(18) 



T w := min {k > N : V(fc) = 0} . 
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Theorem 9 Let the conditions of Theorem 7 be valid, 



max ~Eri} +K < oo, 



E|-0-E G -0r < oo 



and, in addition, 



E]|B| 



< oo for some q > and E \\A\\ < oo. 



Then, as y — > oo 



P (6 > y) ~ C/i/ 



C/ 6 (0,oo). 



(19) 



In particular, 



E<d x < oo 



(20) 



z/ and on??/ i/ x < ft. 

In Section 7 we use Theorem 9 to show that the tail distributions of the 
busy periods of a wide class of branching type polling systems whose input 
parameters and service disciplines vary in a random manner and switch-over 
times are positive, decay at infinity as y~ K for some k > 0. 

4 Auxiliary results 

The proof of Theorem 9 uses Condition T whose validity is not easy to check. 
Wc list here a set of assumptions given in [30] which imply Condition T. 

Let A(A) be the spectral radius of the matrix A. The following statement is 
a refinement of a Kesten theorem from [30] . 

Theorem 10 (see [22]) Let {A n ,n > 0} be a sequence of iid matrices generated 
by a measure Pa with support concentrated on nonnegative matrices and A = 

( a ij)Tj=i = Assume that the following conditions are valid: 

1) there exists q > such that E \\A\\ q < oo; 

2) A has no zero rows a.s.; 

3) the group generated by 

{log A(a„ • • • do) : a n ■ ■ ■ oq > for some n and Oj E supp(J?A)} 

is dense in R; 

4 ) there exists kq > for which 




and 



E \\A\\ Ko log + \\A\\ < oo. 
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Then there exists a n £ (0, k$] such that 

s'(k)= lim ilogE||A l _ 1 .-. J 4 o f = 0. 

n— >oo 77, 

//, in addition, the tuple of m- dimensional vectors {C n ,n > 0} is such that 
the pairs (A n , C„), n e No are iid: (A n , C„) = (A, C) and such that 

P(C = 0)<1, P(C>0) = 1, E||Cf<oo, 

then there exist a constant Kq £ (0, oo) and a continuous strictly positive func- 
tion l(u) on U+ such that relation (14) holds. 

The next lemma provides estimates ol moments of a random walk (see, for 
instance, [29], Theorem 1.5.1): 

Lemma 11 If Xi,i = 1,2,... is a sequence of iid random variables such that 
E \Xi\ p < oo and ~EXi — if p > 1 and M is a stopping time for the sequence 
r„ := X\ + ... + X n then there exists a constant R p £ (0, oo) such that 

E|r M | p < i? p E|X/EAf p / 2vl . 

From now on we agree to denote by K, K\, K2, ... positive constants which 
may be different from formula to formula. 

Lemma 12 ([16]) If a subcritical MBPRE starts by Z(0) = z individuals and 
parameter n in (10) exceeds 1 then for each x £ [1,«0 there exist p = p(x) £ 
(0, 1) and K = K(x) < 00 such that 

E||Z(n)|r<Jrp"||z|f (21) 

for all n £ No- 

The next statement is borrowed from [32], page 351. 
Lemma 13 If conditions of Theorem 6 are valid and, in addition, 

E||B|| 9 <oo and E \\A\\ q < 00 for some q > 0, 
then there exist positive constants K\ and K2 such that for all t > 

P(T > t) < K ie - K2t . 

Let 

Z / (i,j;n,fc) := (Z((i,j; n, k), Z^(i,j;n, k)) 

be the m— dimensional vector of the total progeny alive at moment k > n of the 
j—th. immigrant of type i entering the process at time n, 

m Vi( n ) 
»=i 3=1 
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be the to— dimensional vector of the total progeny alive at moment k > n of all 
the immigrants entering the process at time n, and 



Y(n) :=r?(n)+ £ Z'fafc) 



be the ?n— dimensional vector representing the total progeny steaming from the 
immigrants entering the process at time n. 

Lemma 14 If the conditions of Theorem 7 are valid and, in addition, 

max Er^ +1 < oo 

l<z<m 

then there exists a constant K\ £ (0, oo) such that for all x > 

P(||Y(1)|| >x) <K 1X - K . 

Proof. We have 

oo m ViO-) 

Y(l)=r ? (l) + ^^^Z / (z,j;l,fc + l). 

fe=l i=l j=l 

Hence, observing that, for any fixed i G {1, to} and any j = 1, 2, 77^(1) 

{Z J (i, j; 1, fc + 1), fc = 1, 2, ...} = |Z(fc), fc = 1, 2, ... I Z(0) = e 4 } 
and using Theorem 7 with ^ = 1, i £ {1, to} we see that 



P(||Y(l)||>x) < P(hll>f)+E P E 



i=i \ j=i 



^Z / ( i ,j;l,fc + i; 



fe=l 



> 



2to 



< 



2«E ||r,| 



< 



2 K E||?7 

2 k E||t?| 

x K 



i=l 



E E 

m 

E E 



»?»(!) 

E 

J= i \k=i 



E p E H z ( fc )ll ^ — ^tttI z C°) = 



(i)p £iiz(fc)n 



> 



A 'iE E 



rim 



vfe=l 

(m»h(l))' 



rw?j(l) 



Z(0)=e 4 ;f/ i (l) 



i—1 



The lemma is proved. 
Let 

/i = = min {n : ||V(n)|| > r} 

be the first moment when the size of the population of the MBPIRE exceeds 
level r with the natural agreement that \i = 00 if max„ ||V(n)|| < r. 
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Lemma 15 Under the conditions of Theorem 9 for any e > there exists 
r\ = 7"i (e) > such that 

P ( \\Y(n)\\>ey^<r\ < ey- 

for all r > r±. 

Proof. In view of the estimate 

oo 

^n~ 2 = 7r 2 /6<2 (22) 



n=l 



wc have 



T-l 



P E l|Y(n)||>ey,/i<T 



\n— /x+1 



(OO OO - \ 

E ii Y WII ^{M < « < T} > 6^ 2 £2/ e ^ 
n=l ra=l / 



<EP(HY(n)||I{M<n<T}>!^ 

n=l 
oo 

= EP(M<«<T)P (ll Y(n) || > , (23) 

71=1 

where the last equality is justified by the fact that the event < n < T} is 
specified by V(0), V(l), V(n — 1) while Y(n) depends on the environment 
and immigration after the moment n— 1. Besides, the manipulations with sums 
are correct since the sums run only till T which is finite with probability 1 by 
Lemma 13. 
Clearly, 

P(l|Y(n)ll>^)=P(HY(l)||>^). (24) 



Now using Lemma 14 we have 



P(I|Y(DII>^)< £KyK 



ey \ 2 K Kn 2K 



Thus, 

/ T-l 



P( E HY(n)||> eV>/ i<T] < -^E- 2Kp (M<-<^) 

\n=/*+l / y n=l 

< _^ 2 _ E [t 2k+1 /{ a1 <t}1 < — 

- £ Ky K l y K 

for r > r±(e) since ET 2fc+1 < oo in view of Lemma 13 and /i = fj,(r) t oo in 
probability as r — > oo. 
The lemma is proved. 
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Lemma 16 If 

then for any fixed r 



E(||4||+E|Mir<oo 



E[||V( A1 )|| K /{ / i<T}] <oo. 
Proof. It is easy to check that for any y,w > 0, any e € (0, 1) and x > 

y x < (1 + e)w x + K \y - w\ x , 

where K = K(x,e) := ^1 — (1 + s) l ^ x ^j . Hence introducing a temporary 
notation V[0,n) := (V(0), .., V(n - 1)) we have 

E[||V(n)|r|V(0),..,V(n-l)] 



E [||V(n - l)A n _! + V(n) - V(n - lL^f |V[0, n) 



< (1 + e) || V(n - \\ X + KE [|| V(n) - V(n - 1)A„_ X || x | V[0, n)] 



< (1 + e) || V(n - l)A»-if + 2^E ||77(n)|| 



+2^^ ^E 



Vi(n-l) m 
fe=l i=l 



V[0,n) 



Recalling the definition a,ij(ri) = Ef n (n), set 

m 

and let 



M x (n; i) := E f |ft(n)| !B , M x := max E|/3 l (n)| a: = max E 

Ki<m l<2<m 



Observe that M K < oo in view of condition (16). 
By Lemma 11 for x > 1 we have 



(25) 



(26) 



3=1 



E 



Vi(n-l) m 

fe=i j=i 



V[0,n) 



< 771* £E 



3=1 



Vi(n-l) 

£ fej( n-1 ; fe )- a ij( n - 1 )] 



k=l 



V(n - 1) 



< iW£M x (n " 1, W (»- 1) 



(=i 
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Thus, 



E[||V(n)ir|V[0,n)] < (1 + e) ]|V(n - l^-if + 2 ^ E \\v(n)\\* 

m 

+2 x KR x m x M *( n ~ !> 0^ X/2V1 (" - 1) 

i=l 

= : A K (n- 1). 



(27) 



On the other hand, for x < 1 

E[\\V(n)\\ x \V[0,n)} < E[\\V(n) - V (n)\\ x + \\r)(n)\\ x \V[0,n)} 

< (E[\\V(n) -r,(n)\\\V[0,n)}r + E\\r,f 

< \\V(n-l)A n ^\\ x + E\\ V \\ X <A x (n-l). (28) 

Clearly, 

A K (n-l)7{n-l<M<T} < Q n _i(r) := (1 + e)r K \\A n ^ f + 2 K KE \\ v f 

m 

+2 K KR K m K r K/2yl ^ M K (n - 

i=l 

Observe that by our conditions the expectation E [Q n -i(r)] is independent of 
n and, moreover, 

E[Q„_x(r)] = (l + e)r K E\\A\\ K + 2 K KE\\r,f 



-2 K KR K m K r K/2vl E 



J2M K (n-l,i) 



< (1 + e)r K E \\A\\ K + 2 K KE \\rj\\ K + 2 K KR K m K+1 r K ' 2yl M K := K x < oo. 

Besides, 

||V( M )ir/{M<T} = A K ( M - 1) < Y} < Q„-i(r)- in 



< ^ Q„_ x ( 



A«(m-1) 
||V(n) 



A K (/i-l) 



/i<n<T 



A„(n-1)' 
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Now by (27) and (28) we see that 



E[]|V( M )ir/{M<T}]<E 
l|V(n) 



< E 



A K (n-l) 



/j,<n<T 

/ {T > n} 



IIVWII" 

A„(n-1) 



-i(r) 



/{T>n}E[||V(n)|r|V[0,n)] 



Ve 

^ A K n-l 

n— 1 L v 7 

oo oo 

< E [Q«-i(r)J {T > »}] = ]T E [Qn-i(r)] P (T > n) 

n—1 ri— 1 

oo 

> n) < oo. 

n=l 

The lemma is proved. 
Let 

Z(i,j;l,n) := (Z 1 (i,j;l,ri),...,Z 1 (i,j;l,n)) 

be the to— dimensional vector of the total progeny alive at moment n > I of 
the j—th particle of type i existing in the process at time I. Let, further. 
Z(Z,Z) := V(Z) and 

m Vi(0 

Z(Z,ra) = (Zi(Z,n),...,Z m (Z,n)) := ^ ^ Z(i,j;l,n) 

»=1 3=1 

be the total progeny at time n of all particles existing in the process at moment 
Z < n. Recalling definition (5) we set 



and for n > fj. let 



r(/i )At ) := (V( M ),C M ) = (Z(^ M ),C M ) 



r (M;") : = X! ( Z ( i d;iJ',n),C n ) = (Z(/i,n),C„) 

i=l 3=1 



(29) 



be the inner product of C„ and the vector of the total progeny at time n of the 
V(/x) particles presented in the process at time and 



T-l 



T-l 



r(/x) :=^r(/i ) n) = ^(Z(/i,n),C„). 

tl— fi n—fi 

Denote by C„ (n = 1, 2, . . .) be the a-algebra generated by the tuple 
H (s, A), Hi(s, A), . . . , H n _!(s, A), V(0), . . . , V(n). 



(30) 
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The next lemma shows that for large r the random variable T (/i) = V is, 
in a sense, close to the conditional expectation E [r (//) |C M ] = V(/i)S M (recall 
definition (13)). 

Lemma 17 Under the conditions of Theorem 9 for any e > there exists 
T\ = ?'i (e) smc/i £/ia£ /or all y > 

P (|r ( M ) - (V00, S M )| > £ y; a < T) < [|| V(/i)|r 7 { M < T}] . (31) 

Proof. Evidently, for n > fjL + 1 

(z(/x,n)- v( M )n M ,„,c„) = ^ (z(^,i)n Iin -z(^i-i)n l _i, nj c n ) 

= ^ (z(A*,o-zO*,i- 1)^-1. ni,„c„> 
1=11+1 



which implies 



J2 (Z(fi,n)-V(fi)U^ n ,C T , 



< J2 KzO*,n)-v(M)n^ in ,c r 



E 



£ (z(fx, i) - z(ft, i - i)^_!,n;,„a 



i=JU+l 



OC OC 



< 52\(ZM-Z(l*,l-l)A i -i,U l ,nC n )\ 

l—fi-\-l n—l 



< J2 l|z(M,0-zO*.i-i)^-il 



n ; .„c„ 



£ ||Z( M , Z) - Z( M , ^ - || 
l=tl+l 



(32) 



Here the interchange of summation order is justified since the process having 
V(/i) particles at moment \x < T < oo dies out with probability 1. 
Thus, 

P(|r( / i)-(V( / i),H M )|> £ y;/i<T) 

(OO 
£ ||Z( M , Z) - Z( M , Z - lMi-iH ||H,|| >ey;^<T 
l=ti+l 
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Hence, using (22) we conclude 

P(|r( M )-(V( M ),H M )|> £ y;/i<T|^) 



( C 1 \ 

E I|Z(m, n + /i) - Z( M , n + a* - l)A.+ M -i II l|Hn+ M ll > y E ^ < T | C A 
n—1 n—1 / 

°° / 1 

< \\Z(fi, n + fi)- Z(/x, n + a* - 1)4^-1 II ||H„ +jU || > ^ — ; \i < T ( 



Since Z(/i, n + /1) — Z(/i, n + n — l)A n+p _i and S„ +M arc independent random 
objects on the event /i = k < T < 00 and S n+M = S (see (13)), we get 



£y 1 

P ( ||Z(At,n + A*) - Z(/i,n + At - l^+jt-ill ||H n+M || > y^; M < T 



= jf P (||Z(At, n + n)- Z( M , n + At — e dA ; /i < T|C„) P (||H|| > 

According to Condition T there exists a constant ifi G (0, 00) such that for all 
n > 



_ey_ 

2tn 2 



P ( ||Z( / i,n + At)-Z( / i,n + /i-l)A n+Al _ 1 ||||S„ +A1 || > ^-l ;/i <T 



C, 



< / P(||Z( / i,n + M)-Z(Ax,n + Ai-l)A„ +M _ 1 || Gdt;fi<T\e 



< ^n 2K E [||Z(At, n + ft) - Z(ft, n + At - l)4, +/1 -i f / {At < T} I C M ] . (33) 

Now we consider the cases k < 1 and k > 1 separately. 
For the first case we use for / > At the estimate 



E [||Z(/i, - Z(At, I - 1)4-111* / {M < T} I C„] 



< E 



|Z(Ai,0-Z( / i,Z-l)A i _ 1 || 2 7{Ax<T}|C, 



k/2 



(34) 



Further, we have (recall (26)) 



E 



\\Z(fi, I) - Z( A t, I - l)A,_if 1 {At < T} I C, 



E 



<=i fe=i 3 -=i 



£ E f /3 2 (; - 1)E ^ - 1)/ {At < T} I C M ] 

2=1 

m m 

= £ E f /? 2 (; - 1) own^). < HVOi)^,-! || E Ef/sfcz - 1). 



(35) 
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Thus, for k < 1 and Z = n + /i > /i 

(E [||Z(>, n + m) - Z(/x, n + m - l)4»+ M -i f 1 if* < T l I ^ 



ft/2 



k/2 



< ||v(^)n^ in+ ^_i 

This, in view of the inequality 



s(k/2) = lim (e ||n . 

71— >-00 \ 



,k/2 



^E f ^(n + M-1) • (36) 



,k/2 



1/n 



< 1, 



the first part of condition (17), and relations (35)-(36) leads to the estimate 
P(|r( M )-(V(/i),B |i )|> fi y;i*<T) 



< E 



J2 n 2K \\V(V)\\ K/2 l|n M> „. +M - 1 |r /2 /{ M < T < oo} 



.71=1 



K 



-E 



iiv( M )ir /2 /{ M < t < oo}g^E|in ,„-iir 



/2 



ra=l 



< 7^-72 E [l|V(M)l| K ^{M < T < oo}] < -E [||VGu)fl{/. < T < oo}] 



e y r 



for all r > r (e), proving the lemma for k < 1. 

For the case k > 1 we use Lemma 11 to conclude that for any I > fi 

E [||Z(/i, - Z(/i, Z - l)^_i|| K 7 {/i< T < co} |C„ ] 

m 

< i? K m K M kQ\ i)E [||Z (/i, J - l)f /2vl J < T < oo} |C„ 



< 



i? K m K E [||Z(/i, 



/ - 



ik/2V1 



By Lemma 12 there exist constants K = K (k/2 V 1) <oo and p = p(n/2 V 1) € 
(0, 1) such that 



E 



Z GM - 1)H K/2V1 1 {H < T < oo} \C» < Kp l -»-' ||V(aO|| 



ik/2V1 
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This inequality combined with condition (16) yields for I = n + /i the estimates 
P(|r( M )-(V( M ) 1 H M )|> £ y; M <T) 



KxR K m K - 
< E 



]Tn 2K E 



< 



" Ln=l 

K 1 R K m K+1 M K 



\Z(n,n + (i- 1)|| K/2V1 J2 M «( n + M - 1; i)K 



I {fj, < T < 00} 



E 



™ 2K E [||Z (/i, n + n - 1)|| K/2V1 1 cj I{n < T < 00} 



KiR K m K+1 M K K 

<. Hi 



e K y K 



^n 2K ||V( M )|r/ 2v V- 1 /{M<T<ao} 



£ K y K 



-E 



< £ ^1 /2V1 E til v ^ 11 K 1 < T >] ^ ^ E til ir 1 < T}] 

(the last is valid by selecting r sufficiently large) which justifies the statement 
of the lemma for k > 1. 
The lemma is proved. 



5 The accumulated amount of the final product 

In this section we deduce some estimates related with the total size of the final 
product accumulated in a subcritical MBPIFPRE within a life-period. 
Let 

A y (0) :=e(0)+V(l)J{||r?(l)||>0} 

and let 

m Vi(n) 

A v (n) := ^2 X] <Pi(n;k),n > 1) 
i=i fc=i 

be the total amount of the final product produced by the particles of the n-th 
generation of a MBPIFPRE, 

Aj(n) = A v {n) + ip(n + 1) 

be the total amount of the final product produced by the individuals of the n-th 
generation of a MBPIFPRE plus the final product contributed by the particles 
immigrating at moment n + 1, and 

ro Vi(n+1) / T-1 m Z{ (i,j;n+l,k) \ 

i=l j=l \ fc=n+2 t=l w=l J 

be the total size of the final product contributed to the system by all particles 
immigrating in the process at moment n + 1 and by their progeny. 
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Of interest will be the quantity 

t=i t=i 

i.e., the total amount of the final product produced by the particles of the ra-th 
generation, n > /i, of the MBPIFPRE which belong to the progeny of the V(/x) 
particles of generation fj, and the random variable 

T-l 

Q z ( M ) := J2®z^,n) 

n—fi 

which is equal to the total amount of the final product produced by the V(^t) 
particles of generation fi and their progeny up to the moment of extinction of 
the MBPRE generated by the V(/x) particles of the \i— th generation. 
Clearly, the total amount 

T-l 
n=0 

of the final product accumulated in the system during the life period which 
starts at moment n = 1 admits on the set fi < T the representation 

ft-l T-l 

e = ^A / (n) + ez( M ) + ^r / («)- 

n— n— 

Our aim is to investigate the asymptotic behavior of the tail distribution 
of 6. 

As the first step in solving this problem we compare &z (/•*) with This 
will be done by the arguments similar to those used to demonstrate Lemma 17. 

Lemma 18 Under the conditions of Theorem 9 for any e > there exists 
r = r(e) such that for all y > 

P (|6z (M) - m\ > ey; n < T) < ^-E [||V( M )|| K I {y. < T}] . (37) 

Proof. For ! >fi we have 

m Z t (/i,i) 

e^(M,0-(z(M,0 I c i )=^ (<pt(hk)-E F <p t (i)). 

*=i fe=i 

Now we consider separately the cases k < 1 and k > 1. 
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For k > 1 we use Lemmas 11 and 12 with p = p(k/2 VI) and K = K (k/2 V 1) 
to get 



(\e z (»,i)-{z( f i,i),c l }\> ^ic^ 



A K n 2K 



Z t (M,0 



J2 E M^)-E m (0) 



< (4m) K i? K ^E [||Z( M ,0|| K/2V1 \C,\ |^(Z) - E m (0| K 

< (4m) K i? K ^-% ||V(/i)|r/ 2vl p l ~» E Ef |^(0 - E f ^ t (0| K ■ 

(ey) 



t=i fc=i 

2/,- 



Hence, in view of condition (16) and with I = p + n we have 
P(\B Z (p)-T(p)\>ey;p<T) 
< E 



EPF(|e z (M,M + n)-(Z( M ,/x + n),C /1+n )| > ^|c m )/(a*<T) 

OO 

nv(/i)ir /2vl E" 2 v/(M<T) 



(4m)^ B A-A- l E 



-E 



iv(/z)ir /2vi /( M <T)jE» 2 v i 



< 



N)V f-W2vi) E [ll V Mir j ^ < T )l * ^E[||V( M )f /(„ < T)] 



for all r > r(e). 

To analyze the case k < 1 we apply for I > /j, the inequality 



k/2 



e [|e z (/i, o- (z(/x, 0, Q>r |c„] < (e [|e z (/i, o- (z( M , o, Q}| 2 1 c 

Further, we have 

771 

E[|ez(^,/)-(Z( A/ ,0,Q)| 2 |C A1 ] <™ 2 E E f/3 2 (0E[^(^0 |C„] 

1=1 

m tci 

= m 2 E E f/? 2 (0 (v^n^), < m 2 ||v( M )n^|| E E ^ 2 (0- 

i=l i=l 

Thus, for ft < 1 

k/2 / m \ 

(e [|6zOM)- (z(m,0,c ; )| 2 |c m ]) < m« liv^n^f 72 E E ^ 2 W 



/2 



23 



This combined with the assumption (17) shows that for I = fj, + n 

P(|0 Z (m)-T(m)| >ey;M<T) 
< E 



< -E 



m K Ki 



J2 P (|0^(M, M + »)- (Z(M, M + »), C M+n }| > ^ C M ) / (/i < T) 

oo 

^n 2K ||V(/i)|r /2 ]|n^ +n ir/ 2 J( M <T) 

oo 

||V(/i)|r/ 2 I( M <T)^n 2K E||no, n |r/ 2 



E 



< p^72 E [HV(M)ir J(M < T)] < ^E[||V(/i)ir < T)] 



for all r > r(e). 

The lemma is proved. 

Up to now we have assumed that the initial number of particles and the 
initial amount of the final product are nonrandom. The next two lemmas are 
free of this restriction. 

Let 

JV 



9 (TV) :=£>/(») 



n=0 



Lemma 19 Let a MBPIFPRE be subcritical, and n = 1 be the starting point 
of a life period of the process. If there exists 5 > such that 



E6 K+<> (0) < oo, max Ep? +S < oo, Ei/> K+d < oo, 

l<i<m 

then for any e > there exists r — r(e) such that for all y > yo 

1 



P(e( A1 -l)>ezj; M <T)< 



y 



:+S/2 ' 



(38) 



(39) 



Proof. Let, as earlier, T be the length of the life period in question. Re- 
calling (2) put 



JV 



Q(N,r) :=e(0) + ]T U'(Z + l) + EE^ fc ) 



1=0 



i=i k=i 
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Clearly, 

P(9(iV- l,r) > ey;iV < T) 



< P ( 6 (o) > ^) p (jv < x) + p ^ # + 1) > f J 

(N 7TL r \ 
Z=0 i=l fe=l / 

< f — ^ Ee K+5 (0)P(7V<T) + (7V + l)p(V>> 



ey/ w v ' v y V 3(iV + l) 



(JV + l)rP £>> 



3(JV + l)ry 

< fJL^ F,<d K+5 {0)P(N < T) + (N + 1) ( 3 ^,^ 1) ) E ^ K+5 
In view of Lemma 13 there exist constants -K3 and such that for y > 

E p(^<^)<^- 

Af>if 4 lny ^ 

Now using conditions (38) we have 

oo 

P(6(//-l) >ey; M <T) = E P ( 6 - X ) > £ ^ = # < T) 

JV=1 

< E P(6(iV-l,r) >ey;/i = 7V<T) 

1<N<K 4 Iny 

+ E P (N <T) 

< T~~^+5 E (^V + D^ +1 

l e 2/J l<N<Ktlny 

K 6 ln K+S+2 y 1 

< < 



k+5/2 ' 



The lemma is proved. 



Lemma 20 Under the conditions of Lemma 19 there exists r such that for all 
V>Vo 

P(e>^>T)<-^. (40) 



y 
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Proof. Letting K3 and K4 be the same as in the previous lemma we have 
for y > 2/0 

P(Q>y;Li>T) <P(9(T,r) >y) 

< P (9 {[K A In y] , r) > y) + P (T > K A In y) 

[/f 4 lny] . 

<p(e(o)>|)+ £ pU(/ + i)> 



2/ 



J=0 

[-K4 In y] m r 

2/ 



3 [if 4 In y] 



i=0 i=l k=l 



+ £ EEp(^fe)> 3 mr [4inJ +p(T>X4ln ^ 



T ,ln y 1 
< if 



yK+5 yK + S ' 



The lemma is proved. 



6 Proof of Theorem 9 

Now we are ready to prove the main result of the paper, Theorem 9. First ob- 
serve that by the equivalence of the norms ||-|| and ||-|| 2 and estimates (31), (37), 
(39) and (40), for any e G (0,1/3) one can find r = r(e) such that for all 

V > yo(r,e) 

P (6 > y) < P ((V( M ), S m > > 2/(1 - 3e); ll < T) 



+P (|r(/z) - (V( M ), > ey; » < T) + P (|9 Z ( M ) - T( M )| > ey; M < T) 



+P (9 > y; n > T) + P (9 (p - 1) > ey; Li < T) 



< P «V00, S M ) > y(l - 3e); fx < T) + [||V( M )|| 2 K i < T}] + . 

(41) 

Let l(u) > be the function involved in Condition T. By this condition and 
the independency of V(/i) and S M given [i < T we conclude 



lim sup y K P «V(/z), 3„) > 2/(1 - 3e);^ < T) 

y— >oo 



< lim sup y* 

y—*oo 



s > 



ifo(l - 36)-* 
if (l - 3e)" K E 



P (V(At) G du; /z < T) P 

[ X P(V( M ) edu^<T)\\u\r 2 l(^\-) 

■i||u||,=r VII U ll2/ 

i{^<T} 



(l-3e)2/ \ 

hi 2 y 



i u ii 2 =i 

|V(/i)||^ 



< CO. 
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Since y K P (9 > y) does not depend on r and e, the previous estimate and (41) 
yield 

lim sup y K P (9 > y) < oo (42) 



and, moreover, 

lim sup y K P (9 > y) < K lim inf E 

y— >oo r— »-oo 



|V(a*)J|~I 



v(m) 



/{m<t} 



(43) 



I|v(m)|| 2 

To get a similar estimate from below we use for e > the inequality 

P(e>y) > P(9z( Ai ); M <T) 

> P ((V( M ), H„) > 2/(1 + 2 £ ); M < T) - P (|0 Z ( M ) - r(/x)| > ey;fi < T) 
-P(|rOi)-(V(M),S M )|>ey;M<T) 

Now we select r as large to meet estimates (31), (37) and (39). This gives for 
sufficiently large y > r the inequality 

P (9 > y) > P «V( M ), S„> > y(l + 2e);/x < T) - [|| V(/i)|| 2 K / {/i < T}] . 

Letting 2/ tend to infinity we obtain 



lim inf 2/ K P((V( / x),S A1 )>2/(l + 2e); / i<T) 



(l + 2e)y 



lim inf 2/ K P(V(/ i )Gdu;/ i <T)pf/ F ^ r ,s\ 

J||u|| 2 =r \\IHl2 / ll U 



K (l + 2e)- K E 



I|v(m)II 2 k / 



v(m) 
I|v(m)II 



7{/t<T} 



Let J£T(Z) := inf ue u + K u ) > 0. We know that for sufficiently small e > and an 
appropriate r 



# (l + 2e)- K E 



which implies 



|V(/i)|| 2 K Z 



v(m) 



/{ M <T} 



2 £ E[||V( M )|| 2 /{ M <T}] 



I|V(m)Ii 2 , 

> (A- A-(0(l + 2 e )- K -2e)E[||V( M )|| 2 K /{/ i <T}] >0 



lim inf y K P (9 > y) > 

y-s-oo 



leading in turn to 

lim inf 2/ K P (9 > y) > K lim sup E 



y— >-oo 



l|V(/x)||^ 



Vfr) 
|V(/i)|| 5 



7{ M <T} 



This combined with (43) shows existence of the limit 

V(/i) 



lim E 



l|V(/i)||^ 



l|v(//)|| : 



/{ M <T} 



G (0,00) 
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and gives 



lim y K V (Q>y) = K lim E 

y— too r— >oo 



|v(^)|| 2 K / 



l|V(/i)|| 



I{V<T} 



proving (19). The validity of (20) is now evident. 
The theorem is proved. 



7 Polling systems 



We consider a polling system consisting of a single server and m stations with 
infinite-buffer queues indexed by i G {1, . . . , m} . Initially there are no customers 
in the system and the server waits their arrival at a parking place R. Customers 
arrive to the queues in accordance with a point process whose parameters are 
changing in a random manner each time when the server switches from station 
to station. When a customer arrives to the system, say to station i, the server 
immediately begins the service by visiting the stations in cyclic order (i — > 
i + 1 — m — > 1 — >•••) starting at station i according to a selected service 
policy (to be described later on) and with random positive switch-over times 
between the queues. Later on the initial stage of services (i — > i + 1 — ► • • • — >■ 
m — > • • • — > i — 1) will be called the zero cycle. The subsequent routes of the 
server will be called the first cycle, the second cycle and so on. The server 
immediately moves to its parking place if there are no customers in the system 
to be served. 

Our goal is to investigate various characteristics related with busy periods 
of the system whose work starts by the arrival of a single customer to a station 
i G {1, m} at moment 0. 

To give a rigorous description of the arrival and service processes for the 
system in question we need some notions. Let 



1, ...,m 



be the mixed probability generating functions (m.p.g.f.'s) of (m + 1)— dimen- 
sional vectors (9n, . . . , #i m ; <f>i), where Oij are nonncgative integer- valued random 
variables and <pi is a nonncgative random variable. Denote 



X (s;A) := . . . ,X (m) M)) 



the respective vector-valued m.p.g.f. 
Let, further, 



P W (s;A) 



E 



1 *2 



„Cim p -Ao 



be the m.p.g.f.'s of (m+ 1)— dimensional vectors (Qi, . . . , Q m ; <7j), where Qj are 
nonnegative integer-valued random variables and &i is a nonnegative random 
variable. Denote 

P (s;A) := (pW(s;A),...,p( ro )( S ;A)) 
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the respective vector-valued m.p.g.f. and let H = {(x(s;A), p(s;A))} be the set 
of pairs of the described m.p.g.f. 's. Assume that a probability measure P is 
specified on the natural er-algebra of T-L and let 

(xo(s;A),/9 (s;A)), (xi(s;A), pi{s;\)), (x 2 (s;A),p 2 (s;A)), 

be a sequence of pairs (x«(s;A), p n (s;X)) selected in iid manner from H accord- 
ing to measure P, where Xn{s;X) := n£n^(s;A), . . . , x! m) (s;A)) and p„(s;A) := 

(^Pn \s;X), . . . , ^'"'(s^)^ have the components 



x£ l) M) :=E 



Jn(n) 9 i2 (n) 



, s 6 im (n) e -\<f>i(n) 



^ } (s;A) 



E 



c C>i(n) &,(n) 
*1 *2 



>(™)g- A(T i(™) 



In the present paper we investigate such polling systems whose arrival and 
service procedures of customers meet the following two conditions. 

Branching property. When the server arrives to station i for the ro-th time 
and find there, say, ki customers labelled 1, 2, hi , then, during the course of 
the server's visit, the arrival of customers is arranged in such a way that after 
the end of each stage of service of customer j (the number of such stages may be 
more than one if the service discipline admits feedback) the queues at the system 
will be increased by a random population of customers (9n(n,j), . . . , Oi m (n,j)) , 
where 6u(n,j) is the number of customers added to the Z-th. station, and, in 
addition, a final product of size (j>i(n,j) > will be added to the system. It is 
assumed that the vectors (0u(n,j),...,9i m (n,j);(j)i(n,j)), j = 1,2, ...,ki are 
iid and such that 



x! l) (s;A). 



Note that Xn\&',X) is a random m.p.g.f. and, therefore, the parameters of 
the polling system are changed from cycle to cycle in a random manner. 

Immigration property. Upon departing from the i— th station for the 
n— th time the server switches to the queue existing at station (i + 1) (mod to). 
The process of switching takes a random time at the end of which the queues at 
the system will be increased by a random population of customers (Qi (n) , . . . , Q m 
where Cu( n ) is the number of customers added to the /-th station, and, in addi- 
tion, a final product of size Ui (n) > will be added to the system. It is assumed 
that the vectors (Cii(^)i • ■ • 5 C»m( n )j &i{ rl ')) > n = 1, 2, ... are iid and such that 



, s <n(n) s C i2 (n) 



In the sequel we call the described polling systems as the branching type 
polling systems with final product and random environment (BTPSFPRE). 

Polling systems are rather popular subject of investigations in queueing the- 
ory (see surveys [25] , [35] and [43] for definitions and more details) . For instance, 
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polling systems possessing the branching property in which the probability gen- 
erating functions hn\s), n = 1,2,... arc nonrandom and identical for all 
cycles were considered in [12], [36], [41] and quite recently in [4] and [5]. These 
models cover many classical service policies, including the exhaustive, gated, 
binomial-gated and their feedback modifications. 

Polling systems with input parameters and service disciplines changing in a 
random manner and (or) depending on the states of systems were investigated 
by the fluid method in [9]- [11]) and by the method based on the construction of 
appropriate Lyapunov functions in [13] - [15]. The present paper may be viewed 
as a compliment to article [16] in which the BTPSFPRE with zero switchover 
times were analyzed by means of MBPRE without immigration. 

Before we pass to general results, consider two examples of BTPSFPRE. 

Let 7+ = {T} be the set of all probability distributions of nonnegative random 
variables, 

T™ xl™ :={(7i,...,T m ) x {I 1 ,...,I m ):T i eT + ,I j eI +> i,j = l,...,m} 

be the set of all 2m-dimensional tuples of such distributions, M e = {£} be the 
set of all m x m matrices £ = (£ij)™ = i with nonnegative elements, and M 7 = 
{r} be the set of all m x (m + 1) matrices F = (7y )™ 1:)= g with nonnegative 
elements such that 

m 

y^7y = l, i = l,...,m, max 7,0 > 0. 

L — ' l<2<m 
3=0 

Let P be a measure on the Borel a— algebra of the space 

S := Me x Me x M-y x T™ x Tff. 

Example 21 (motivated by [15]). Consider a polling system with m stations 
and a single server performing cyclic service of the customers at the stations. 
Assume that initially there are no customers in the system and the server is 
located at a parking place R. Assume that given the idle system the flow of cus- 
tomers arriving to station I is Poisson with, say, a deterministic rate £; . When 
the first customer appears in the system the server selects a random element 

(£ ,£o J ,r ,To xl ) es with 

£0 = (e y (0))™ =1 ,^=(4(0))"; =1 , ^ = (^(0))™^, 

To = (Tio, T m o) , Io = (Iio, Imo) 

and immediately starts its zero service cycle (i — > i + 1 (modm) 
i — l(modm)) adopting the gated server policy with random switch-over times. 
Namely, the server serves all the customers that were queueing, say at sta- 
tion i when the server arrived and then proceeds to the next (in cyclic or- 
der) station with switch-over times Cj(0) distributed according to Iio(x) := 
P(c7;(0) < x). During the switch-over period <Ti(0) new customers arrive to 
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the system according to independent Poisson flows with intensities given by the 
vector (efi(O), ef 2 (0): ■•■j £ fm(0)) (some of the components may be equal to zero). 

For the period while the server performs the batch of services at station i 
new customers arrive to the system according to independent Poisson flows with 
intensities given by the vector (£ii(0), £^(0), ...,£j m (0)) (some of the components 
may be equal to zero) and the service times of customers are iid and distributed 
according to Tio(x) := P (Ti(0) < x). Each served customer either goes to station 
j 6 {1, . . . , m} with probability 7y(0) or leaves the system with probability 7io(0) 
independently of other events. Besides, after the end of each service stage of a 
customer the customer contributes to the system the service time at this stage 
as the final product. 

It is assumed that given (£q, £q , To, To x Io) the service times, switch-over 
times and the arrival process of new customers are independent. 

The subsequent routes n = 1, 2, . . . have the same probabilistic structure spec- 
ified by the tuples 



f 

T„ 



( e ij( n ))Tj=l ' £ n= ( £ ij( n )Y" 



(7ij(»))2=i t ,-= o; 



• . 1 1 n 



with only difference that at the beginning of cycle n > 1 there is a possibility to 
have more than one customer in the system. 

One can show (see [16]) that this system possesses a branching property in 
which the m.p.g.f. x„(s;A) has the components (in our setting and the service 
time of customers as the final product) 



X«M) : =E 



Jn(n) 9 a (n) 
*1 b 2 



i(™)g- AT i(™) 



A + £ £ij .(n)(l- Sj ) I ,(44) 



where 



t m (A) 



e~ Xx dT in (x), 



and an immigration property with p„(s;A) having components 



p^M) 



E 



,Cii(") Ci2(n) 



. s Ci m (n) e ~\cn(n) 



= ^°°exp j- ^A + fj4(n)(l- Sj )j zjdffeOc). (45) 

Example 22 (compare with [15]). Consider the same polling system as earlier 
but assume now that at each station the server adopts the exhaustive server 
policy: it serves all the customers that were queueing at the station when the 
server arrived together with all subsequent arrivals up until the queue at this 
station becomes empty and then switches over to the next station. 
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Let 



w irl (A) := 



1 - lu{n) 
1 - lii{n)t ln {\) 



tin (A) 



(46) 



and 



yin[s) 



7io(n) + E 7 ^7ijH 



1 - 7*; («) 

It is know (see, for instance, [16]) that the system above possesses a branch- 
ing property with c/)„(s;A) whose components are unique solutions of the equa- 
tions 



xl l) (s;A) = y in (s)w in A + £ e«(n)(l -*,•) + e«(«)(l - 

The immigration property for this system fulfills with p„(s;A) specified by (45). 

Now we come back to the general case. We assume that the server may start 
the service of a customer of any type arrived to an idle system immediately 
and would like to study the distribution of the total size of the final product 
accumulated in the system during a busy period of the server. For this reason 
the law of arrival of customers to an idle system plays no role for the subsequent 
arguments. The only assumption wc need that the probability of arrival two or 
more customer to an idle system simultaneously is zero. 

Set 

h®(a) := x W(s;0),i=l,2 > ...,m 
and for n = 0,1,2,... introduce m.p.g.f.'s 



F«(s;A) = E 



J,l(n) £,2(«) 
4>j t> 2 



. q (im(n) -\ipi(n) 



and p.g.f.'s 



/^(s) = E 



„&i(n) Ji2(n) _ _ _ f im (n) 



by the equalities F^fs; A) := Xn ( s jA), 

FW(s; A) := X » • • • , 'uFg+V (s; A) , . . . , (s; A) ; a) , * < m, (47) 
and/^s) :=/il TO) (s), 

/W(s) := /#> («!,..., Si, (s) , . . . , /("»> (s)) , i < m. 



(48) 



Further, set 

iW(s) :=pW(s;0) J i = l I 2,...,m 
and for n = 0, 1, 2, . . . introduce m.p.g.f.'s 

G„(s;A) := E [s? {n) sf {n) ■ ■ ■ j^We^W 



/ft 

11^ (*!.•••> s - F i i+1) (b; A) , .. . , ^i m) (s; A) ; a) (49) 
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and p.g.f.'s 



g n (s) := E \s r ( 



(s),...,/< m) (s)) 




(50) 



We would like to describe conditions on the branching type polling system 
under which power moments of the amount of the final product accumulated in 
the system during its busy period are finite or infinite. Note that letting the 
final product 4>i(n,j) be the service time of the j-th customer served during 
the n-th visit of the server to station i and <7i(n) be the switch-over time from 
station i to station (i + 1) (modm) after the n-th visit of the server to station 
i we provide conditions under which the tail distribution of the length of the 
busy period of a polling system decays, as y — > oo, like y~ K for some k > 0. 

Our results are based on an important statement revealing connections be- 
tween the behavior of certain characteristics of the busy periods of BTPSFPRE 
and the related characteristics of life-length periods of a MBPIFPRE. To for- 
mulate this statement we introduce the notion of generalized busy period as 
follows. 

If a busy period of a BTPSFPRE starts by the arrival of a single customer 
at station J then, after the end of the busy period generated by this customer 
we continue to follow the performance of the system (with the described laws 
of arrival of customers within switch-over times, the service disciplines and ac- 
cumulating of final product) up to the first visit of the server to station J when 
there are no customers in the system. 

Definition 23 The time interval between the start of service of the first cus- 
tomer and the moment when, for the first time, the system is idle after the end 
of the switch-over time from station J — 1 to station J is called a generalized 
busy period. 

Denote by Mj the length of the respective generalized busy period. 

Since a generalized busy period includes a possibility to have an idle system 
at the moment when the server arrives to station J — l(modm) but have cus- 
tomers at the moment when the server arrives to station J (which arrive to the 
system during the switch-over time J — 1 — > J), the length of a generalized busy 
period is not less than the length of the standard busy period. 

In the sequel we assume without loss of generality that J = 1. 

Theorem 24 If a generalized busy period of a BTPSFPRE starts by the arrival 
of a single customer at station 1 then the distribution of the amount of the 
final product accumulated in the system to the end of the generalized period 
coincides with the distribution of the total amount of the final product produced 
in a MBPIFPRE within a life period which starts at moment by the birth of a 
single particle of type 1 and where the joint distribution of the number of direct 
descendants and immigrants and the amount of the final product produced by 



33 



particles of different types of the k-th generation is given by the vector-valued 
m.p.g.f. 's 

H fc (s;A) := (s;A), . . . , *f l) (s;A); G fc+1 (s;A)) , k = 1, 2, . . . , n 

specified by (47) and (49). 

To prove this theorem one should modify in a natural and evident way the 
proof of Theorem 3 in [36] and we omit the respective arguments. 

The MBPIFPRE described in Theorem 24 will be called the associated for 
the BTPSFPRE. 

Now we apply the obtained earlier results for MBPIFPRE to our polling 
system. 

Let H n := (hij (n)) i - =:l be the matrix with elements 



rf(s) 



ds i 



and A n := (dij{n)) i =1 be the matrix with elements 



oy(n) := 



5/^ } (s) 



9s, 



= E f „^(n). 



(51) 



Then in view of (48) a mj -(n) = hmj(n),j — 1,2, ... ,m and for z < m 

m 

aij(n) = hij(n)I {j < i} + h ik (n)a kj (n). 

fe=i+l 



(52) 



For i = 1, . . . , m introduce auxiliary matrices 
••• 



/ 1 




ffW := 



■■• 






V o 



hi(i-i)(n) hu(n) h l{%+1) {n) 













him (jl) 

o 




1 



where for each i the elements of the matrix H n are located only in row i of the 
matrix H„. It is not difficult to check by (52) that 



A - H {1) H {2) ■ ■ -H (m) 



(53) 
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Further, let C n := (Ci(n), . . . , C TO (n)) be a random vector with components 
dF®(l;\) 



Ci{n) := 



dX 



A=0 



E Pn ¥>i(n), i = 1, . . . ,m 



and let c„ := (ci(n), . . . , c m (n))' be a random vector with components 



Ci(n) 



dA 



A=0 



E^ n 0i(n), i = 1, . . - ,m. 



Then, by (47) C m (n) = c m (n) and, for i < m 

m 

Ci(n) = Ci(n)+ ^ h ik (n)C k (n). 
fc=t+i 

Hence we get C„ = (E 1 — 1 c„ where 

^:=(^-(n)/(i<j))™ =1 

is the upper triangular matrix generated by H n . 

Let, further, L n := (Zjj ("))™ =1 be a random matrix with elements 



= Ez„Cij(n), 



and let 



Pi(n) 

Then, by (49) and (50) 

dg n (s 



dp£ J (l;A) 



cZA 



Bjin) 



A=0 



E p „(Ti(n), i = l,...,m. 



m m 



and 



A, = E G „^(n) = 



dG„(l,A) 



rZA 



A=0 



= y^Pi(n) + hk{n)C k (i 



\i=l k=i+l / 

The next two statements are easy consequences of Theorem 9 . 

Theorem 25 Assume that the MBPIFPRE associated with a BTPSFPRE is 
subcritical and satisfies conditions of Theorem 9. Then there exist constants 
C\,C2 € (0, oo) such that 



C lV - K < P (Q P > y) < C 2 y~ K , y > y Q . 



(54) 



In particular, if the final product of any customer is its service time then the 
tail distribution of the length Op of a busy period of the system satisfies (54)- 
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Corollary 26 Under the conditions of Theorem 25 EOp < oo if and only if 

x e (0, «). 

Proof of Theorem 25. We associate with the initial polling system an aux- 
iliary polling system obtained from the initial one by excluding certain customers 
and the final product contributed by them. We split the exclusion process into 
several stages. 

At the first stage we exclude from the initial system all the customers (and 
the final product contributed by them) which arrive to the system during the 
switch-over times of the server. We call these customers the customers of the 
first level. 

At the second stage we exclude from the initial system all the customers (and 
the final product contributed by them) which arrive to the system during the 
service times of the first level customers. We call these customers the customers 
of the second level. 

At the k-th stage we exclude from the initial system all the customers (and 
the final product contributed by them) which arrive to the system during the 
service times of the (k — l)-th level customers. And so on. 

Denote by $ - the total amount of the final product accumulated in the 
auxiliary branching type polling system. It is clear that Op > <!>. As shown in 
[16], under the conditions of Theorem 25 P (<& > y) ~ Ky~ K as y — > oo. Hence 

P (6p > y) > P ($ > y) > C l2 T K , y > y . 

This proves the needed bound from below. 

To get the desired estimate from above note that Op < Q ex t where Q e xt is 
the total amount of the final product accumulated in the analyzed branching 

type polling system within the generalized busy period. By Theorem 24 <d ex t = 
© where © is the total amount of the final product produced in the associated 
MBPIFPRE within a life period which starts at moment by the birth of a 
type 1 particle. Using now Theorem 9 we obtain 

P (Op > y) < P (0 > y) ~ C iy ~ K , y y oo, 

which gives the desired estimate from above. 

For completeness we formulate the following statement concerning polling 
systems with a > 0. 

Theorem 27 Assume that the MBPIFPRE associated with a BTPSFPRE is 
such that its underlying MBPRE satisfies conditions of Theorem 5 with a > 
and, in addition, 

P [ min Em > ) > 0. 

Vl<i<m j 

J/0p is the total size of the final product accumulated in the BTPSFPRE during 
a busy period then P (Op = oo) > 0. In particular, if the service time of any 
customer at any station is positive with probability 1 then the busy period of the 
BTPSFPRE is infinite with positive probability. 
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Proof. As shown in [16] under the conditions of Theorem 27 P ($ = oo) > 
where 4> is the same as in the proof of Theorem 25. Hence the desired statement 
follows. 
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